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The two-photon exchange contribution to the single spin asymmetries with the spin orientation 
normal to the reaction plane is discussed for elastic electron-proton scattering in the equivalent 
photon approximation. In this case, hadronic part of the two-photon exchange amplitude describes 
real Compton scattering (RCS). We show that in the case of the beam normal spin asymmetry, 
this approximation selects only the photon helicity flip amplitudes of RCS. At low energies, we 
make use of unitarity and estimate the contribution of the nN multipoles to the photon helicity 
flip amplitudes. In the Regge regime, QRCS approximation allows for a contribution from two pion 
exchange, and we provide an estimate of such contributions. We furthermore discuss the possibility 
of the quark and gluon GPD’s contributions in the QRCS kinematics. 
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I. INTRODUCTION 

Recently, the new polarization transfer data for the 
electromagnetic form factors ratio Ge/Gm raised a 
lot of interest for the two photon exchange (TPE) physics 
in elastic electron proton scattering. These new data ap¬ 
peared to be incompatible with the Rosenbluth data pj. 
A possible way to reconcile the two data sets was pro¬ 
posed [|j, which consists in a more precise account on 
the TPE amplitude, the real part of which enters the ra¬ 
diative corrections to the cross section (Rosenbluth) and 
the polarization cross section ratio in a different manner. 
At present, only the IR divergent part of the two pho¬ 
ton exchange contribution, corresponding to one of the 
exchanged photon beeing soft, is taken into the experi¬ 
mental analysis Q. Two model calculations exist for the 
real part of the TPE amplitude and they quali¬ 

tatively confirm that the exchange of two hard photons 
may be responsible for this discrepancy. In order to ex¬ 
tract the electric form factor in the model independent 
way, one has thus to study the general case of Compton 
scattering with two spacclike photons. These two photon 
contributions are important for the electroweak sector, as 
well. 

In view of this interest, parity-conserving single spin 
asymmetries in elastic ep-scattering with the spin orien¬ 
tation normal to the reaction plane regain an attention. 
These observables are directly related to the imaginary 
part of the TPE amplitude and have been an object of 
theoretical studies in the 1960’s and 70’s Q. By ana- 
lyticity, the real part of the TPE amplitude is given by 
a dispersion integral over its imaginary part. Therefore, 
a good understanding of this class of observables is ab¬ 
solutely necessary. Recently, first measurements of the 
beam normal spin asymmetry B n * 1 have been performed 
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1 In the literature, also the A n notation for beam normal spin 
asymmetry or vector analyzing power was adopted. 


in different kinematics [gj. 

Though small (several to tens ppm), this asymmetry 
can be measured with a precision of fractions of ppm. 
Before implementing different models for the real part, 
where an additional uncertainty comes from the disper¬ 
sion integral over the imaginary part, one should check 
the level of understanding of the imaginary part of TPE. 
These checks have been done for the existing data. Inclu¬ 
sion of the elastic (nucleon) intermediate state only [9J led 
to negative asymmetry of several ppm in the kinematics 
of SAMPLE experiment but was not enough to describe 
the data. The description of the beam normal spin asym¬ 
metry within a phenomenological model which uses the 
full set of the single pion electroproduction m did not 
give satisfactory description at any of the available kine¬ 
matics. Especially intriguing appears the situation with 
the SAMPLE data with electron lab energy Ei a b = 200 
MeV, which is just above the pion production thresh¬ 
old where the theoretical input is well understood. On 
the other hand, an EFT calculation without dynamical 
pions iia was somewhat surprisingly very successful in 
describing this kinematics for B n . This success suggests 
that to the given order of chiral perturbation theory, the 
role of the dynamical pions for this observable might be 
not too large. Finally, a recent calculation within the 
leading logarithm approximation appeared, where only 
few dominant multipoles were used as input (a. Sur- 
prizingly, the authors of E2 were able to describe all the 
experimental data at very different beam energies and 
scattering angles, while the full calculation of [l0| failed 
to describe any of them. 

Though even at low energies the situation with the 
imaginary part of TPE amplitude is by far not clear, 
an attention has to be paid to high energies, as well, 
since the dispersion integral which would give us its real 
part, should be performed over the full energy range. 
Due to relative ease of measuring B n within the frame¬ 
work of parity violating electron scattering, new data 
from running and upcoming experiments yd} will stim¬ 
ulate further theoretical investigations of this new ob¬ 
servable. A calculation of B n in hard kinematics regime 
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at high energy and momentum transfer was performed 
recently in the framework of generalized parton distribu¬ 
tions (GPD’s) and resulted in asymmetries of ~ 1.5 ppm. 

H; 

Since a ppm effect measurement at high momentum 
transfers is an extremely complicated task, the forward 
kinematics seems more favorable. For this kinematics, a 
calculation exists E3 , where an observation is made that 
the contribution of the situation where the exchanged 
photons are nearly real and overtake the external electron 
kinematics is enhanced as ln 2 (— t/m 2 ) ~ 100, with m the 
electron mass and t < 0 the elastic momentum transfer. 
Making use of the optical theorem, the authors obtained 
estimates of B n in this kinematics as large as 25-35 ppm. 
The calculation of (Tb] appeared afterwards showed that 
this result is not adequate, since the squared log term can 
only contribute with the Compton scattering amplitudes 
with the photon helicity flip. Indeed, this conclusion was 
confirmed in Ref. G2 , the corrected version of M . it 
was noticed that the In 2 (—t/m 2 ) term should vanish in 
the forward kinematics due to gauge invariance, and the 
leading term is governed by a single (though still large) 
log term, In (—t/m 2 ) ~ 10. The predictions for B n shifted 
correspondingly to ~ —(4 — 6) ppm and agree with the 
preliminary data from HAPPEX experiment 11 .‘fl j. 

In order to provide an estimate of B ni we use the equiv¬ 
alent photon or quasi-real Compton scattering approxi¬ 
mation which is caused by the hard collinear kinematics, 
responsible for this ln 2 (— t/m 2 ) enhancement. In this 
approximation, the leading contribution comes from the 
kinematical region with both exchanged photons are al¬ 
most real. The hadronic tensor is taken at this kinemati¬ 
cal point and can be taken out of the integration. For the 
hadronic tensor, we adopt the most general real Compton 
scattering amplitude and demostrate that the contribu¬ 
tion of the cross section (i.e., photon helicity conserving 
amplitudes) vanishes in the QRCS approximation. The 
remaining contributions are related to the photon helic¬ 
ity flip amplitudes. We provide the calculation of B n at 
low energies, where nN intermediate states are expected 
to be the dominant contributions. In this kinematics, the 
asymmetry can be related to the pion photoproduction 
multipoles, and we discuss the relative contributions of 
different multipoles. At higher energies, we discuss for¬ 
ward kinematics, that is Regge regime, and note that 
the combinations of the RCS amplitudes appearing in 
the expression of B n are related to 27 t exchange in the 
f-channel. We provide an estimate of such a contribu¬ 
tion but conclude that these are negligibly small. We 
furthermore discuss the hard regime and discuss possible 
contributions which might be enhanced by the large log¬ 
arithms originating from the QRCS kinematics. 

The article is organized as follows: in SectionlTTl we de¬ 
fine the kinematics, general ep-scattering amplitude and 
the observables of interest; in Section mu the two pho¬ 
ton exchange mechanism and the photons kinematics is 
studied; the equivalent photons or quasi real Compton 
scattering approximation and its implementation for the 


case of B n is given in Section m we present our results 
in Section El and conclude with a short summary. 


II. ELASTIC ep-SCATTERING AMPLITUDE 


In this work, we consider elastic electron-proton scat¬ 
tering process e{k) + p(p) —> e(k') + p(p') for which we 
define: 


p = 

p + p 

2 


K = 

k + k' 


2 


q = 

k — k! = p' — p, 

(1) 


and choose the invariants t = q 2 < 0 2 and v = (. P -K)/M 
as the independent variables. M denotes the nucleon 
mass. They are related to the Mandelstam variables s = 
(p + k ) 2 and u = (p — k') 2 through s — u = AMv and 
s + u + t = 2 M 2 . For convenience, we also introduce 
the usual polarization parameter e of the virtual photon, 
which can be related to the invariants v and t (beglecting 
the electron mass m): 


v 2 — M 2 t(1 + r) 

v 2 + M 2 t( 1 + r) ’ 1 ’ 

with r = — t/(4M 2 ). Elastic scattering of two spin 
1/2 particles is described by six independent amplitudes. 
Three of them do not flip the electron helicity [3|, 


1 no flip 


—^{k’^^k) (3) 

U (P) (GmY - F 2 — + U(p ), 


while the 
thus have 
0 : 


other three are electron helicity flipping and 
in general the order of the electron mass m 


_ m e~ 
fUp ~ M -t 




u(k')u(k) ■ u(jp') ( F a + F 5 — ) u(p) 


+ F 6 u(k')^u{k) ■ u{p')j 5 u(p) 


( 4 ) 


In the one-photon exchange (Born) approximation, two 
of the six amplitudes match with the electromagnetic 
form factors, 

Gfr n M = c M {t\ 

F 2 Born (u,t) = F 2 (t), 

= 0 (5) 


2 In elastic ep-scattering, the usual notation for the momentum 
transfer is Q 2 = — q 2 but we prefer the more general notation t 
to avoid confusion with the incoming and outgoing photon vir¬ 
tualities in forward doubly virtual Compton scattering we will 
be concerned in the following. 
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where Gjf(t) and F 2 (t) are the magnetic and Pauli form 
factors, respectively. For further convenience we define 
also Ge = Gm — (1 + t)F 2 and F\ = Gm — F 2 which 
in the Born approximation reduce to electric form factor 
Ge and Dirac form factor F \, respectively. For a beam 
polarized normal to the scattering plane, one can define 
a single spin asymmetry, 


B 


n 


Cp — Cj 
f7 T + CT| ’ 


( 6 ) 


FIG. 1: Two-plioton exchange diagram. 

The imaginary part of the two-photon exchange (TPE) 
graph in FiglU is given by 


= e 2 


|fci | 2 rf|A;i |<iO fcl 


2£i(27t) 3 * 


MTi'ffl + m )lv u 


QlQl 


W^(w\Qi,Q 2 2 ) 


( 10 ) 


where a j (oq) denotes the cross sesction for an unpolar¬ 
ized target and for an electron beam spin parallel (an¬ 
tiparallel) to the normal polarization vector defined as 


cm _ 



[fcxfc'A 
\kxk'\) ’ 


(7) 


normalized to (S • S) = —1. Similarly, one defines the 
target normal spin asymmetry A n . It has been shown 
in the early 1970’s Q that such asymmetries are directly 
related to the imaginary part of the T-matrix. Since 
the electromagnetic form factors and the one-photon ex¬ 
change amplitude are purely real, B n obtains its finite 
contribution to leading order in the electromagnetic con¬ 
stant a em from an interference between the Born ampli¬ 
tude and the imaginary part of the two-photon exchange 
amplitude. In terms of the amplitudes of Eqs. (I3I4II . the 
beam normal spin asymmetry is given by: 


— _ £ )' v ^^ + T { t ^M + £ Ge) 

■ |rGMlnrF 3 + G_ElmF 4 + Fi— ImF 5 |. (8) 


For completeness, we also give here the expression of 
target normal spin asymmetry T n 3 in terms of invariant 
amplitudes: 


T n = {tG" m + eG%) 

+ t) F-[ lmF‘2 — FJmFi 
2e „ „ \ v 




+ 


— G B -G M -LnF 3 


(9) 


III. TWO PHOTON EXCHANGE 



3 Also An notation for target normal spin asymmetry exists in the 

literature. 


where W^ u {w 2 , Qf, Q 2 ) is the imaginary part of doubly 
virtual Compton scattering tensor. Q\ and Q\ denote 
the virtualities of the exchanged photons in the TPE di¬ 
agram, and w is the invariant mass of the intermediate 
hadronic system. We next study the kinematics of the 
exchanged photons. Neglecting the small electron mass 
and using the c.m. frame of the electron and proton, one 
has: 


Ql, 2 = 2|fc||fci|(l - COS01 >2 ), (11) 

with |fc| = s 2 ^ = k the three momentum of the incom¬ 
ing (and outgoing) eletron, |fci| = \J{ S ~ W 2 A ~ ) 2 ~ 

that of the intermediate electron, and cos 02 = 

cos 0 cos 0i + sin 0 sin 0i cos <j>. The kinematically al¬ 
lowed values of the virtualities of the exchanged photons 
(the restriction is due to the fact that the intermediate 
electron is on-shell) are represented by the internal area 
of the ellypses shown in Fig. [3 The ellypses are drawn 

Accessible values of Qj 2 , Q,, 2 



FIG. 2: Kinematically allowed values of the photon vitualities 

Qi, 2. 


inside a square whose side is defined through the external 
kinematics (k) and the invariant mass of the intermedi¬ 
ate hadronic state ( w 2 or k\), while the form solely by 
the scattering angle. If choosing higher values of the 

mass of the hadronic system w 2 < s, it leads to scaling 

_ 2 

the size of the ellypse by a factor of 5 _ . In the limit 
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w 2 = (y/s — m e ) 2 , the ellypses shrink to a point at the 
origin and both photons are nearly real. This is not a 
soft photon (IR) singularity, however, since the real pho¬ 
tons’ energy remains large enough in order to provide the 
transition from nucleon with the mass M to the interme¬ 
diate state X with the mass w. Instead, the intermediate 
electron is soft, /cf « (m e , 0 ), therefore this kind of kine¬ 
matics does not lead to an IR divergency which can only 
occur if the intermediate hadronic state is the nucleon 
itself. In the following we are going to study this kine- 
matical situation in more detail. 


IV. QUASI-RCS APPROXIMATION 

We consider the kinematical factors under the integral 
over the electron phase space of Eq. m at the upper 
limit of the integration over the invariant mass of the 
intermediate hadronic state, w —> yfs — m e : 

hi. I ^ I ^ J__L ( 12 ) 

E 1 Q\Ql Ak 2 E\ mAh 2 ' V ’ 

In this range of the hadronic mass w, the exchanged pho¬ 
tons are real, and the contribution of real Compton scat¬ 
tering (RCS) to the 27 -exchange graph is enhanced by a 
factor of 1/m (it is not a singularity since B n has a factor 
of m in front). 

We next rewrite the contraction of the hadronic and 
the leptonic tensors in Eq. m identically as 

l^W^(W 2 ,QlQ 2 2 ) = 0,0) 

+l$W tv (8,0,0) 

+w [W^(W 2 , Ql Q\) - W^(s, 0,0)] , (13) 

where we expanded the leptonic tensor l llv = 
with 

1% = mu'w^u 

(14) 

In Ea. (ll3ll . only the first term does not vanish in the 
QRCS limit, while the second and third are at least linear 
in k\. In the limit W 2 —> s, one has f,\ —+ my 0 , which is to 
be compared to the term m -1 in 1° . Though on the first 
glance on the first glance they are of the same order, the 
additional 7 -matrix picks an extra electron mass to from 
one of the projectors jt 4 - to when performing the sum¬ 
mation over electron spins. Quasi-real Compton scatter¬ 
ing (QRCS) approximation consists in assuming the first 
term to be dominant due to the kinematical enhancement 
under the integral and in neglecting the second one. In 
general, this kind of contributions coming from QRCS 
kinematics will allways be present in the full calculation, 
since the second term in Eq. m is constructed in such 
a way that the resulting integral is regular at the QRCS 
point. In the following, the QRCS approximation will be 
used. Hence, the hadronic tensor can be taken out of the 


integration over the electron phase space. The remaining 
integral is 


In = 


d 3 ki 


2E 1 (2n) 3 Q 2 Q 2 


2 1 


The result of the integration reads: 


(15) 


h 


1 

-32n 2 t 





+ 8 Sp 



,( 16 ) 


where Sp(x) is the Spence or dilog function, Sp(x ) = 

- /q 1 y ln(l - xt) with Sp( 1) = 7 r 2 / 6 . In the high energy 
limit, —> 1 , we recover the result of Ref.jTlj. For 
the details, we address the reader to the Appendix. This 
In 2 (^|) factor serves as the justification of the QRCS 
approximation. The full result may be represented as 
a “Taylor expansion” in powers of ln(^|) ~ 10. By 
studying the QRCS approximation we show in a model 
independent way that the leading term in this expansion 
is quadratic. Since this leading term is large, we can 
(under ceretain kinematical conditions) neglect further 
terms. 

The RCS tensor may be taken for instance in the basis 
of Prange 0 or, equivalently of Berg and Lindner [l9l| . 

w rcs = + ^ 2 ) + + mo 

P^n v - n^P’ v 

+-- nmr 

plfj, v I n !iptu 3 

+ pX," N - < 17 > 


with the vectors defined as P' = P — = 

p u K a qp such that P'-K = P'-n = n- K = 0. 
The amplitudes Bi are functions of v and t. This form 
of Compton tensor is convenient due to the simple form 
of the tensors appearing in Eq. EJ. 

Before contracting the leptonic and hadronic tensors, 
we notice that the amplitude Bj can only contribute to 
the invariant amplitude F 6 , since it contains the N'j 5 N 
structure. Fg does not contribute at leading order in 
to to neither obserevable of interest, therefore Bi will 
be neglected in the following. The remaining tensors in 
Eq. (Q Til are symmetric in indices p,v. 

ImM^CS = e 2 mIo u'uW^ cs g^- (18) 

Finally, we identify different terms in Eq. Cl with the 
structures, together with amplitudes Fi t ...e parametriz¬ 
ing elastic ep-scattering amplitude in Eqs. cm and 
find for the invariant amplitudes for the elastic electron- 
proton scattering in the QRCS approximation: 

Ini G^ RCS = Tm p Q RCS = Im pQ RCS = 0 (19) 
ImF 4 Qi?CS = -MtI 0 lm(Bi + B 3 ) (20) 

Tm P 9 RCS = -M 2 tI 0 Im(B 2 +B 4 ) (21) 
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Hr, 


We obtain for B n in the QRCS approximation: 

B n = mtlQyJ2e{l — £)Vl + t {jG\[ + eG^) 

■ {G E lm(Bi + B 3 ) + vF 1 lm(B 2 + B 4 )} • (22) 


The result of Eq. (El is obtained by only using the 
assumption that the QRCS kinematics dominate the in¬ 
tegral in Eq. m- The combinations of the RCS am¬ 
plitudes appearing in the final result can be further¬ 
more expressed in terms of the helicity amplitudes of 
real Compton scattering. With these latter defined as 

Tx'X'M N , - 4T4 w* cs > one has 20® 


B\ + B 3 = — 


1 


V~t 
2 M 




M 4 - su 12 ’ 12 
2 M 1 


(23) 


2 M i / _ \ 

B2 + B4. = -—7 —-— (T -i 1 -1 1 T-1 1 -11 ) 

2 s + M 2 


+ 


M 4 — su s — M 2 1 3- _1 5 


(24) 


ImT 


0 


— 1—- 1 - 
-L 0 , J. 0 


= sin — 167rqg. v / sy^(fc + l) 2 


fe >0 


X [|R/c +| 2 — |R ( fc +1) _| 2 ] 


0 


x F(—k, k + 2, 2, sin 2 —) 


0 


= - sin 3 — 87r^ ^ 

fc>i 

x [|Rfc+| 2 — |R ( fc +1) _| 2 ] 


(25) 

k 2 (k + l) 2 (fc + 2) 2 


12 


0, 


x F(—k + 1, A: + 3,4, sin 2 —) 


(26) 


ImTji _ x i = sin 2 cos -^- 8 nqi T \/~s 


k(k + l) 2 (fc + 2) 


x Re 


fc>i 


-®fc+^4fc+ — 5( fc+1 )_R( fc+1 )_ 

. 2 e> 


x F(—fc + 1, fc + 3,3, sin —), 


(27) 


where q n is the c.m. pion three-momentum and the hy¬ 
pergeometric function is defined as 

1 1 , abx , a(a + 1)6(6 + 1) z 2 , ocA 

na.b.c.x) = I + - 1 T+ c(c + 1) 2!+.(28) 


As can be seen, only photon helicity-flipping ampli¬ 
tudes enter the final result for B n in the QRCS approxi¬ 
mation. This is the main result of this work. 


We keep only few first multipoles in this infinite se¬ 
ries, namely Ao+, A 1+ , B 1 +, A 2 _, B 2 _ which obtain 
their leading contributions from threshold pion produc¬ 
tion, A(1232) and Di3(1520) resonances. The result for 
B n reads: 


V. RESULTS 


In this section, we consider the impact of the QRCS 
contributions onto the beam normal spin asymmetry 
in different kinematics: low energies ( ttN intermediate 
states), high energies and forward angles, i.e. Regge 
regime (27 t exchange in the f-channel) and hard regime 
(handbag diagrams and two gluon exchange). 


A. Low energies: 7 tN multipoles 


B n = SirmqT, 


st 


\J 2er (1 — e)\J 1 + t 


(s — M 2 ) 2 tG 2 m + eG 2 e 


- [|^ o +| 2 - l ^ r - l 2 + 4(|^ 1+ | 2 - |^ 2 _| 2 )] 
6 v / s ; 


M 

3 


-FiR e(B 1 + A* 1 + -B 2 _A* 2 _) 


0 


— sm — 
2 2 


- f ■) 


(29) 


(4(|A 1+ | 2 -|A 2 _| 2 ) + |R 1+ | 2 -|R 2 _| 2 ) 


+2 ( — F 2 - —F\ 
\M M 


■Re(B 1 + A* 1 + ^B 2 _A*_)] 


Above the pion production threshold, the imaginary 
part of the RCS helicity amplitudes can be related to the 
pion photo- and electro-production multipoles. These 
relations read mm for the three amplitudes entering 


The helicity multipoles are related to the electromag¬ 
netic multipoles E t+ ^ + M i+ as 

Bq+ = Aq+j Mi- = Ai_, (30) 
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and for l > 1 
Ei+ 

M l+ 

E (i+ 1 )- 


l + l 

1 

ITT 

i 


Ai+ + t;Bi+ 


A 1 + 2 B 
A i+ -x— £>i+ 


Mi 


(i+i)— 


l + l 

1 

7+T 


1 + 2 


Hi+i)- 


B, 


(i+i)— 


A(i+ 1)- + 2' B ( i + 1 )- 


(31) 


It is informative to consider the contributions of Eq+ 
and Mi + which are dominant at low energies. Neglecting 
other multipoles, we get: 


B n = — 87 rmq. 
Elab 


St 


\J2e{\ - e)\/l + t 

1 rGl + sGl 


M 


(s - M 2 ) 2 

F 2 -Fi^j [\E 0+ \ 2 + \M 1+ \ 2 ] 

0 (y/s ~ M ) 2 \y/s + M 


(M Fl+3sin : 2 


2yfsM 


M 


Fo — F-i 


■\M 1+ \ 2 


(32) 


The first term in the brakets is negative for Ei a b < 
Et ~ 0.45 GeV for the proton target, and allways nega¬ 
tive for the neutron target. Due to the overall minus sign, 
the beam normal spin asymmetry necessarily obtains a 
positive contribution from the threshold pion production. 
Moving towards the A(1232) resonance position, we see 
that the dominant contribution now comes from the sec¬ 
ond term. The factor multiplying the |Mi + | 2 term is 
allways positive for the proton and negative for the neu¬ 
tron. Therefore one expects that the A resonance give a 
large negative contribution to B n on the proton target, 
and a positive one on the neutron target. We use the 
MAID2003 multipoles as input to Eo . 112911 and present 
in Fig. 0 the energy dependence of the beam normal 
spin asymmetry for the proton target over the resonance 
region. This result can be compared to the results of 
Refs. [H[T3. If confronted to the full calculation .10|, 
we find agreement with their findings in the QRCS ap¬ 
proximation. In Ref. Pj . the authors do nominally the 
same approximation as we do, keeping the ln 2 (Q 2 /m 2 ) 
and the RCS part of the hadronic tensor only but ar¬ 
rive to the same sign of the contributions of the Eq + and 
M\ + multipoles for the proton. They are furthermore 
able to describe both low energy backward angle data 
from SAMPLE and intermediate energy and angle data 
from MAMI within the same approximation. It has been 
shown in 0 that the QRCS approximation does work 
well at backward angles (in the sense that it represents 
the dominant part of the full integration range) but drops 
short at forward angles, where the exchange of at least 
one hard virtual photon is important (for very forward 
angles, see 01. Therefore, even if the QRCS approxima¬ 
tion did work in the forward regime, one should not take 




1.2 


1.4 


1.2 


1.4 


W„ N (GeV) 


W„ N (GeV) 


FIG. 3: Beam normal spin asymmetry for the reaction e+p —> 
e+p in the QRCS approximation at fixed c.m. scattering an¬ 
gle 0cm = 120° as function of the 7 tN invariant mass W n N ■ 
The contribution of the Eo+ nrultipole (upper left panel), E\+ 
and Mi+ (upper right panel), E 2 - and M 2 - (lower left panel) 
are shown separately. The full result within the QRCS ap¬ 
proximation is shown in the lower right panel. 


this success too seriously, since the neglected double or 
single virtual Compton scattering effects are important. 
We show the angular distributions of B n on the proton 
target for different values of the electron lab energy in 
Fig. a In Fig. 0 we display the energy distributions of 
B n on the neutron over the resonance region. 

It is necessary to note that the resonance region does 
not seem to be favorable for the QRCS approximation. 
In this approximation, the value of the hadronic am¬ 
plitude is taken at the largest energy available and is 
not integrated over the full spectrum. Therefore, if one 
goes above a resonance position, it is only tail of the 
resonance that contribute, which causes the asymmetry 
to drop quite fast (see the sharp resonance behaviour 
at A(1232), Di3(1520) in Fig0. In a full calculation, 
also intermediate energies of the hadronic system do con¬ 
tribute and the asymmetry does not have such a sharp en¬ 
ergy dependence. So, it is preferable to use the QRCS ap¬ 
proximation in the region where the energy dependence 
is rather smooth. Furthermore, the quality of the QRCS 
approximation is function of the scattering angle, as well. 
It works better at backward angles and worse at forward 
angles. The reason for this is clear: taking the limit of 
the soft intermediate electron, fci « 0 and neglecting the 
dependence of the hadronic and leptonic tensors on k\ 
corresponds to performing Taylor expansion of the inte- 
gration result in powers of In and neglecting all the 

2 O 2 

terms beyond the first leading log term ~ In . The 
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©cm ( de g) ©cm ( de g) 


FIG. 4: Beam normal spin asymmetry on the proton as func¬ 
tion of c.m. scattering angle at different values of lab energy: 
200 MeV (upper left panel), 300 MeV (upper right panel), 
570 MeV (lower left panel) and 855 MeV (lower right panel). 
The different lines represent Eo+ multipole contribution (thin 
solid lines), E\+ and contribution (dashed lines), E 2 - 
and M 2 - contribution (dotted lines), and the sum of all (thick 
solid lines) 


larger Q 2 , the better the approximation works. This be¬ 
haviour was first observed in Ref. 0 - Instead, results 
of Ref. [2| suggest that the loop integral in B n is com¬ 
pletely dominated by the QRCS region in all the kine¬ 
matics (backward regime for SAMPLE data and rather 
forward regime for MAMI data). 


B. Regge regime: 27 t exchange 

The calculation of Ref. 0 estimates B n as 

B n ~ cr t ot ln(-f/m 2 ), (33) 

which is based on taking the RCS amplitude in the ex¬ 
act forward limit. If one allows the momentum trans¬ 
fer to be non-zero, one however obtains the contribution 
from photon helicity flip amplitudes, which is the result 
of this work. Though it is logical to expect that these 
amplitudes be suppressed at low values of t, there is a 
competing factor of ~ In ^, which can be of the order 
of several hundreds. In this subsection, we provide an es¬ 
timate of such a contribution at forward kinematics. The 
combinations of the RCS amplitudes B 1 + B 3 , B 2 + B 4 in 
Regge regime are known to have the quantum numbers 
of a scalar isoscalar exchange in the t-channel which was 
successfully described by cr-meson 20j or equivalently, 
by two pion exchange in the t-channel [21 ■ Since the 


FIG. 5: Beam normal spin asymmetry for the reaction e + 
n —> e + n in the QRCS approximation as function of the 7 tN 
invariant mass W-kn- Notation as in Fig. [3 


effective cr-meson exchange does not result in a non-zero 
imaginary part in the s-channel, one should use the two 
pion exchange mechanism accompanied by multiparticle 
intermediate state in the s-channel. In this section we 
estimate B\ + B 3 and B 2 + B 4 through ttN and pN inter¬ 
mediate states contributions within the “minimal” Regge 
model for n and p photoproduction j25j where reggeized 
description of high energy pion production is obtained by 
adding the t-channel meson exchange amplitude and (in 
the case of yp reaction) s-channel Born amplitude which 
is necessary to ensure gauge invariance. The reggeization 
procedure naturally leads to the substitution of each f- 
channel Feynman propagator by its Regge counterpart, 
T^T ->• with 


_ -ML _:_ ( 34 ) 

sin770^(1) T(l + a 7r (t)) ’ 

with an(t) = a4(t — m 2 ) and c>4 = 0.7 GeV -2 . Gauge 
invariance requires the s-channel piece to be reggeized in 
the same way, i.e., to be multiplied by (t—?n 2 )P£ J (a w (t)). 
For compactness, we will use the shorthand I20j . 




B\ + B 2 + v(B 2 + B 4 ) 


(35) 


for the combination of the Bi ’s which enters the result for 
the asymmetry. Here, we list the results of the calculation 





















of Ai: 


Im Al N = 


Im A^ n 


O m 2f2 

K 71 -(RJ-M 2 AJ), 


M(s - M 2 ) 


(36) 




M 2 


2 M{s-M 2 ) [s-M 2 
s — 3 M 2 + 5?n 2 

+ - 2 - 


rCff 


+ M 2 


K 


where C n = 2\/2 Me ^ N with f 2 nn /4tt = 0.08, and 
C p = 2V2 with / pw7 = 0.103 [H. Where 

possible, m n and Q 2 were neglected as compared to s, M, 
m p in order to simplify the final expression. The scalar 
integrals Aq, Bq, Cq for both nN and pN contributions 
are given in the Appendix. 

The result of Eq. (pi leads to negligibly small contri¬ 
bution to the B n , of the order of 10 -2 ppm for energies 
between 6 and 45 GeV, which is to be compared to ~ 5 
ppm from Ref. E3- There are three suppression factors 
which are responsible for this small result. Firstly, the 
Reggeization procedure which leads to suppression in en¬ 
ergy v a '” W. Secondly, the amplitude A± is defined as the 
combination -^[B 1 + B 3 + v(B 2 + -B 4 )], and the singu¬ 
larity j in the definition only cancels if taking the gauge 
invariant combination as described in l25l . Therefore, we 
obtain an additional suppression factor of t. The other 

feature of the result of Eqs. is that both contri- 

2 m 2 

butions are suppressed by factors A - 2 and S _^ I2 , re- 
spectivelyly. In the case of the pion, it is interesting to 
observe this fact in view of somewhat surprising success 
of the effective field description of the SAMPLE data 
point on B n without pion contribution. This might give 
a hint that the pion contribution to B n at low energies 
is suppressed by the pion mass, if calculated to the same 
order. 


C. Hard regime: GPD’s 

Recently, beam normal spin asymmetry was studied in 
the hard scattering regime 0 in the framework of the 
handbag mechanism. In this picture, the process factor¬ 
izes into the hard part (electron scattering off an on-shell 
quark with the exchange of two spacelike photons) and 
the soft part which is parametrized by the generalized 
parton distributions (GPD’s). However, one cannot ac¬ 
cess the QRCS kinematics here, since it is impossible for 
an on-shell quark to absorb (emit) a hard real photon and 
remain on-shell. As it has been noticed before, the QRCS 
kinematics selects the RCS amplitudes with the photon 
helicity flip, i.e. the helicity changes by two units. In the 
context of the GPD’s, helicity amplitudes with the helic¬ 
ity changed by 2 units were considered in several works. 
In the collinear kinematics, this requires that scattering 
occurs on a parton of spin 1, that is gluon. The cor¬ 
responding gluon helicity flip GPD’s were introduced in 


- / ^e iXx <p'S'\F^(~n)F^(~n)\pS> 

X J Z7T Z Z 

p([/x a]- z^/3) 

= H Tg (x, 0 N'ip', S') - Q — - N(p, S ) 

+Et 9 (x,0N'(p',S') - q M f 2 N(p,S) (37) 


where Ht b (x,£), £t s (x,£) are the double helicity flip 
gluon GPD’s, q = p' — p, and the matrix element is 
taken between the initial and final nucleon states with 
the momenta p and p', and spins S and S', respectively. 
The integration is performed along the light-cone vector 
= ( 1 , 0 , 0 , 1 ). x and £ are the longitudinal parton mo¬ 
mentum fraction and skewness parameter, respectively. 
Furthermore, \pct\ and [v(3\ are antisymmetric pairs of 
indices, while (...) means symmetrization and removal 
of the trace. These operations are essential since the 
product of operators must transform as irreducible rep¬ 
resentations of the Lorentz group. This requirement rules 
out the possibility to see these GPD’s in the QRCS kine¬ 
matics, since it is exactly the trace of the RCS tensor 
that contributes to B n in the QRCS regime, c.f. Eq. G3- 
So, our conclusion is that in the QRCS kinematics, con¬ 
tributions from generalized parton distributions are ruled 
out. 


VI. SUMMARY 

In summary, we presented a calculation of the beam 
normal spin asymmetry. This observable obtains its lead¬ 
ing contribution from the imaginary part of the two pho¬ 
ton exchange graph times the Born amplitude and is di¬ 
rectly related to the imaginary part of doubly virtual 
Compton scattering The resulting loop integral obtains 
a large contribution from the kinematics when both ex¬ 
changed photons are nearly real and collinear to the ex¬ 
ternal electrons. We adopt the QRCS or equivalent pho¬ 
tons approximation which allows to take the hadronic 
part out of the integral and to perform the integra¬ 
tion over the electron phase space analytically. For the 
hadronic part, we use the full real Compton scattering 
amplitude and show, that only photon helicity flipping 
amplitudes contribute in this observable in the QRCS 
approximation. At low energies, we relate this helicity 
flipping Compton amplitude to the irN multipoles and 
discuss the contributions of different multipoles. 

At high energies and forward scattering angles (Regge 
regime), we provide an explicite calculation which is due 
to two pion exchange in the f-channel. The resulting val¬ 
ues of the asymmetry are of the order 10~ 2 ppm for the 
energies in the range 6 — 45 GeV. We conclude that the 
double logarithmic enhancement does not dominate B n 
in forward regime since it comes with helicity-flip Comp¬ 
ton amplitude which highly suppresses this behaviour. 

Finally, we consider possible contributions to B n in the 
hard kinematics, among which contributions from quark 
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and gluon GPD’s and show that such contributions are 
ruled out in the QRCS kinematics. 
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VII. APPENDIX A: INTEGRALS OVER 
ELECTRON PHASE SPACE 


In this section we present calculation of the integrals 
over electron phase space appearing in the QRCS approx¬ 
imation. First we calculate the scalar integral /q, 


h 



k\dk\ 

2Ui(2tt 3 ) 


dklk 1 

(k- fci) 2 (fc' - /fci) 2 ’ 


(38) 


where the upper integration limit kthr corresponds to 
the inelastic threshold (i.e. pion production), kthr = 


_ m 2 


4s . We next introduce integration 

over the Feynman parameter A. = J 4 ^ a+ ^l a y x y . We 
chose the polar axis such as k\ ■ (k—xq) = k±\k—xq\ cos 0i 
with |fc — xq \ 2 = k 2 + x(l — x)t and perform angular inte¬ 
gration. We furthermore change integration over dk\ to 
integration over dimensionsless 2 = E\jE 


do 




dz 


m 2 4m 2 (-1 \2 

t V 1 z ) 


(39) 



4^(1-^ +y^ 2 -^ 

^(1 - z )* - y / Z 2 - 


To perform the integration over the electron energy, 
we follow here the main details of the calculation in Ap¬ 
pendix A of Ref. 0- The result reads 


Io 


1 

— 327 T 2 t 



f _ t B thr \ 

\ TO 2 E 2 J 


+ 8 Sp 



,(40) 


where Sp(x) is the Spence or dilog function, Sp(x ) = 
- £ y l n (l — x t) with Sp( 1) = 7t 2 /6. In the high energy 
limit, —> 1, we recover the result of Ref.p~5j. 


VIII. APPENDIX B: SCALAR INTEGRALS 
FOR HELICITY FLIP AMPLITUDE 


The vector and tensor integrals can be reduced to the 
scalar ones by means of standard methods |27|. The re¬ 
maining integrals to be calculated are the two, three, and 


four-point scalar integrals. Here we are only interested 
in the imaginary part of these, therefore there are only 
three integrals with non-zero imaginary part: the two- 
point integral 


Cl = Im 


d 4 pr- 


(27r) 4 p 2 — to 2 (P + K — p n ) 2 — M 2 


87T J~i 


(41) 


with \p n \ = — rrO 2 ; the three-point one 


B 


7T 

0 


f d 4 pn 1 _ 1 _ 

7 (27t) 4 pI - m 2 (k - Pit ) 2 - ml 
1 

’ (P + K -p n ) 2 -M 2 

1 . 2 E„ 

87r(s — M 2 ) m v ’ 


(42) 


and finally, the four-point integral: 

d 4 p-x 1 


Al = Im 


1 


(27r) 4 (k-pt r) 2 - m 2 p 2 - m 2 
1 1 
' (V - p „) 2 - ml {P + K- pv? - M 2 

1 


8 nQ 2 (s — M 2 + m 2 ) 
1 


1 + 


4 mlE 2 
Q 2 Pl 


: In 


1 + 


4 m^E 2 
Q 2 pl 


+ 1 


, 4 mlE 2 1 

1 + Q2 P 2 1 


(43) 


These integrals should however be reggeised as de¬ 
scribed in Section El by substituting the Regge propaga¬ 
tor instead of the Feynman one. Denoting t\ = ( k—p n ) 2 
and Q = ( k’ — p^) 2 the momentum transferred by the 
pions in the f-channel, we have for the reggeized version 
of scalar integrals: 


0 ci) R = J dn^h-mDv? K(u)) 

■(t 2 - (44) 


( B o) R = ^2^ J dflnih-ml) (45) 

P*K(<i))^K(Q)), 


(Al) R = ^ ^Jdn n V?Mh))V?Mt 2 )) (46) 


Similarly, in the case of the p-exchange in the s- 
channel, the integrals with non-zero imaginary part are: 


C? 


T f d 4 p p 1 _1_ 

J (27t) 4 p 2 —m 2 (P + K — p p ) 2 — M 2 

1 \Pp\ 

87r -yfs ’ 
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B p 

n o 



d 4 p p 1_1_ 

(2tt) 4 p 2 p - m 2 p (k - p p ) 2 - ml 
1 


(P + K - pp) 2 - M 2 
1 . 2 E p (s - M 2 ) 

87r(s — M 2 ) Mm 2 p 


with \p p \ = yj (s AI 4 ^ m ——- m 2 , and a 2 = E 2 m 2 p - 
EE p (m 2 — ml) + \{m? p — ml) 2 . If neglecting the pion 

Mm 2 . . . 

mass, a = 2s ' . I he reggeization procedure is the same 
as for 7 tN- intermediate state. 




= Irn 


d 4 p p 


(27r) 4 {k - pp) 2 - m 2 p 2 - m 2 


1 


1 


(fc' - Pp) 2 - m 2 (P + K — pp) 2 - M 2 
1 1 


8ttQ 2 (s -M 2 + m 2 ) ^Ji + _t 

• In 


4cr 2 


1 + ^ + 1 

1 + 


(47) 
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